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Motivated by recent Fermi surface and transport measurements on LaNiOs, we study the Mott 
Metal-Insulator transitions of perovskite nickelates, with the chemical formula RNiOs, where R is a 
rare-earth ion. We introduce and study a minimal two-band model, which takes into account only 
the Qg bands. In the weak to intermediate correlation limit, a Hartree-Fock analysis predicts charge 
and spin order consistent with experiments on R=Pr, Nd, driven by Fermi surface nesting. It also 
produces an interesting semi-metallic electronic state in the model when an ideal cubic structure 
is assumed. We also study the model in the strong interaction limit, and find that the charge and 
magnetic order observed in experiment exist only in the presence of very large Hund's coupling, 
suggesting that additional physics is required to explain the properties of the more insulating nick- 
elates, R=Eu,Lu,Y. Next, we extend our analysis to slabs of finite thickness. In ultra-thin slabs, 
quantum confinement effects substantially change the nesting properties and the magnetic ordering 
of the bulk, driving the material to exhibit highly anisotropic transport properties. However, pure 
confinement alone does not significantly enhance insulating behavior. Based on these results, we 
discuss the importance of various physical effects, and propose some experiments. 
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I. INTRODUCTION 

The Mott Metal-Insulator Transition (MIT) is a cen- 
tral subject in the physics of correlated electron phenom- 
ena and transition metal oxides.^ The perovskite nicke- 
lates, RNiOa, where R is a rare earth atom, constitute 
one of the canonical families of materials exhibiting such 
an MIT. One of the most interesting features of the nick- 
elates is the charge and spin ordering in the insulating 
state, which is relatively complex yet in the ground state 
is robust across the entire family. ^^ The explanation of 
this ordering is still in many ways controversial. While 
the MIT in bulk nickelates is an old subject, the topic 
has been reinvigorated recently by attempts to grow thin 
film heterostructures and observe unique quantum con- 
finement effect s.^^^ In this paper, we revisit the problem 
of the MIT and ordering in the nickelates, both in bulk 
and in heterostructures, from a very simple theoretical 
viewpoint. 

We begin by summarizing some salient features of the 
nickelates. First, as the rare-earth ionic radius decreases, 
the MIT temperature increases. Starting from R=La 
which is metallic at all temperature, R=Pr, Nd have fi- 
nite MIT temperatures Tmit = 120K, 180K respectively 
(R=Eu has the highest MIT temperature Tmit = 480K) 
and finally R=Lu is insulating at all temperatures. This 
trend is understood due to the increasing distortions in- 
troduced in the smaller rare earth materials, which in- 
crease the Ni-O-Ni bond angle and hence reduce the 
bandwidth. In the materials R=La, Pr, Nd, the elec- 
trons can therefore be understood as more itinerant and 
bandlike, while they are increasingly "Mottlike" for the 
smaller rare earths. 

Second, at low temperature, all the nickelates display a 
magnetic ordering pattern with an "up-up-down-down" 



spin configuration which quadruples the unit cell rela- 
tive to the ideal cubic structure. ^^ This pattern coex- 
ists with a "rock salt" type charge order - what is ac- 
tually observed is expansion or contraction of the oxy- 
gen octahedra - which alternates between cubic sites. 
Such charge order must indeed always be present for this 
magnetic state, on symmetry grounds, and can there- 
fore be considered to this extent as a secondary order 
parameter. ^^ Interestingly, for the more metallic nicke- 
lates both charge and spin order appear simultaneously, 
consistent with this view, while for the more insulating 
nickelates, R=Eu, Ho, the charge ordering occurs in- 
dependently in an intermediate temperature insulating 
phase without magnetism. ^'^'^^ 

A variety of microscopic physical mechanisms have 
been proposed for the nickelates. A naive view of the 
material would be to consider the nickel d electrons only, 
occupying the nominal Ni^+ valence state which would 
place one electron in the e^ doublet, which is degenerate 
with cubic symmetry. Early studies attributed the com- 
plex spin pattern to orbital ordering, perhaps induced by 
Jahn- Teller or orthorhombic distortions that split the e^ 
degeneracy. However, no Jahn-Teller distortion was ob- 
served, and it was later suggested that orbital degeneracy 
is removed by a separation of charge into Ni^+ and Ni^+ 
states (an extreme view of the charge order), which have 
no orbital degeneracy. ^^ This was attributed to strong 
Hund's rule exchange on the nickel ion,^^ but phonons 
may also be involved. However, the observed and ro- 
bust magnetic ordering is not so natural in this picture. 
Another question mark is raised by spectroscopic mea- 
surements, which seem to observe a significant of Ni^+ 
occupation, suggesting that a model with holes on the 
oxygens may be more appropriate.^^ In this paper, we 
reconsider the mechanisms for spin and charge ordering 



in the nickelates, and specifically highlight the distinc- 
tions between an itinerant and localized picture of the 
electrons. Our main conclusion is that, at least for the 
R=La, Pr, Nd materials where the MIT transition tem- 
perature is low or zero, and in which a broad metallic 
regime is observed, the itinerant picture is more appro- 
priate. We summarize the main content of the paper 
below. 

The contrasts between the aforementioned models are 
really sharp only deep in the Mott limit, in which or- 
bital degeneracy, ionic charge, and Hund's rule versus 
superexchange are clearly defined and distinct. In an 
itinerant picture, the precise atomic content of the bands 
is not in itself important, but rather the physics should be 
constituted from a model of the dispersion of the states 
near the Fermi energy and the interactions amongst these 
same states. In this view, the observed ordering may be 
considered as spin and charge density waves (SDWs and 
CDWs), and are tied to the Fermi surface structure. Re- 
cent soft X-ray photoemission^^ indeed observed large 
fiat regions of Fermi surface in LaNiOs, which appear fa- 
vorable for a nesting-based spin density wave instability. 

Specifically, in Sec. II we introduce a minimal two- 
band model for the electronic states near the Fermi en- 
ergy in the nickelates. While it is easiest to motivate 
such a model from the naive view of Ni^+ valence states 
- which is questionable, as noted above - it can be 
considered just as the simplest phenomenological tight 
binding Hamiltonian which can produce electronic bands 
with the appropriate symmetry, in agreement with LDA 
calculations^^. Within this model, the crucial parame- 
ter controlling the shape of the bands is the ratio of the 
second neighbor to first-neighbor d-d hopping. With a 
small and reasonable ratio, the large closed Fermi sur- 
face observed in experiments and LDA calculations is 
reproduced. ^^'^^ In addition, the same fermiology reason- 
ably explains the resistivity. Hall effect and thermopower 
measurements on LaNiOs ^^, as well as the main features 
of the optical conductivity below 2eV}^ We study the 
effect of interactions in this model by a simple random 
phase approximation (RPA) criterion for the spin density 
wave instability, and by more detailed Hartree-Fock cal- 
culations, in Sec. III. These mean-field type approaches 
are, we believe, reasonably appropriate for the itinerant 
limit. Interestingly, we find that the same hopping ra- 
tio which reproduces the experimental Fermi surface also 
turns achieves nearly optimal nesting, which further sup- 
ports the itinerant view. The Hartree-Fock calculations 
then predict the phase diagram as a function of spin- 
independent and spin-dependent interactions, which we 
include microscopically by Hubbard U and Hund's rule 
Jh couplings in the tight-binding model. 

We find that the Hartree-Fock calculations produce 
two possible explanations for the observed spin and 
charge ordering in the more itinerant nickelates. The- 
oretically, these two scenarios can be best understood by 
considering a hypothetical ideal cubic sample (the real 
materials undergoing MITs are orthorhombic even in the 



metallic state). In such a sample, we obtain two distinct 
insulating ground states, characterized by "site centered" 
and "bond centered" SDWs. If it occured within an oth- 
erwise cubic sample, the bond centered SDW would have 
equal magnitude of moments on all sites, and would not 
induce charge ordering. In the site centered SDW, charge 
ordering is present, and there would be a vanishing mo- 
ment on one rock salt sublattice. In real orthorhom- 
bic samples, the bond centered SDW will be driven off- 
center, and charge order is induced. The latter off-center 
SDW appears most consistent with experiment. It is also 
the most favorable SDW state in the Hartree-Fock calcu- 
lations, and dominates in the regime of relatively small 
Jh coupling. 

For completeness, in Sec. IV we study the two-band 
model in the strong coupling limit, in which U and/or 
Jh are much larger than the bandwidth. In this limit, 
we find that an insulating state with charge order con- 
sistent with experiment can be obtained, but only for 
very large Hund's exchange, Jh/U > 4. The magnetic 
order is found to be either ferromagnetic or of the site 
centered SDW type. While the latter is quite close to 
what is observed in experiment, it does not appear fully 
consistent, and moreover the requirement of such large 
Jh to stabilize a charge ordered state seems to reaffirm 
the unphysical nature of this limit. 

After this detour to strong coupling, we return to 
the reasonably successful model and Hartree-Fock ap- 
proach, and apply it to finite thickness slabs in Sec. V. 
This provides a minimal and highly idealized model for 
a nickelate film. We find that quantum confinement 
leads to substantial changes of the nesting properties of 
ultra-thin slabs. The predicted consequences are modi- 
fied magnetic ordering compared to bulk and and highly 
anisotropic transport properties. One result we do not 
find from this calculation is a substantial enhancement 
of the Mott insulating state in films of just a few mono- 
layers, a phenomena for which there is gathering ex- 
perimental evidence. ^^'^^^^'^^ We take this as evidence 
that the putative Mott insulating state in ultrathin LNO 
films is driven not only by confinement but by additional 
interface-sensitive effects. 

Finally, we conclude in Sec. VI with a discussion of 
experiments, models, and some open issues. In particular 
we discuss the role of oxygen 2p orbit als, and a possible 
physical mechanism behind the insulating state. We also 
describe some experimental probes of the Mott transition 
which may help to distinguish different mechanisms. 



II. TWO-BAND MODEL AND NESTING 
PROPERTIES 

The simplest tight-binding model for the nickelates is 
constructed based on the naive Ni^+ valence. In this ionic 
configuration, the only partially occupied orbitals are the 
two members of the e^ doublet, containing one electron. 
We consider the hopping through the neighboring oxygen 
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FIG. 1: Fermi surfaces for the tight-binding model. In (a) and 
(b), we show the conduction and valence band Fermi surfaces, 
respectively, for i' It — 0. For larger t' jt^ the conduction band 
Fermi surfaces become large and hole-like, as shown in (c) and 
(d) ioi t' /t = 0.15. The approximate nesting in the latter case 
is indicated schematically in (d). 



p states (cr-bonding) as dominant, and treat it as virtual. 
This leads to strongly directional hopping, described as 
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where z,j are site indices, a, 6 = 1,2 are orbital indices 
for 2z^ — x^ — y'^ and x^ — y'^ respectively and a =^, \, 
is the spin index. Comparison with LDA band calcula- 
tion and with the experimentally measured Fermi sur- 
face indicates that the nearest-neighbor hopping t and 
the next-nearest-neighbor hopping t' with a-type bond- 
ing is the most dominant. In detail, t^^^n = t(j)^^(j)^^ and 

K^±^±o = t\rA + <^') whe^e ^- = i-h "#) 

2' 2 / 



= (1,0) are the orbital wavefunctions for 



the 2x^ — y'^ — z^, 2y'^ —x^ — z^ and 2z^ —x'^—y'^ a-bonding 
orbitals along the three axes. We estimated t' /t ~ 0.05 
by fitting our tight-binding model with LDA band cal- 
culation, while the best fits to experimentally measured 
Fermi surface gives t' /t ~ 0.15.^^'^^ 

The range 0.05 < t' /t < 0.2 reasonably explains 
the observation of an hole-like Hall coefficient but an 
electron-like thermopower in LaNiOs.^^'^^ This appar- 
ently contradictory behavior of the Hall conductivity and 
thermopower arises from the mixed electron and hole 
character on the e^ Fermi surface. The model also pre- 
dicts inter-band optical spectral weight in reasonable cor- 
respondence with experiment at low energy (less than 
2eV).ii 

We now examine the Fermi surface in more detail in 
search of nesting tendencies. Fig.l shows representative 
Fermi surfaces obtained from the tight-binding model as 
a function of the ratio t' /t. With increasing t' /t from to 
0.15, the topology of the large Fermi surface is changing 
as seen in Fig. 1(a), (c) and (d). In the absence of next- 
nearest-neighbor hopping t' /t = 0, the conduction band 
Fermi surface has an open topology as seen in Fig.l (a). 



With increasing t' /t^ this Fermi surface becomes closed, 
comprising a large "pocket" centered at the zone corner. 
In the intermediate range (especially 0.1 < t' /t < 0.2), 
the pocket resembles a cube, as seen in Fig. 1(c) and (d) 
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FIG. 2: Zero frequency spin susceptibility for the tight- 
binding Hamiltonian for t' /t — 0.05,0.1,0.15, as a function 
of momentum k in the cubic Brillouin zone. Note that for the 
best nested situation, t' /t — 0.15, the susceptibility is sharply 
peaked close to the wavevector 27r(|, |, |). 



((d) shows both valence band and conduction band Fermi 
surfaces). Contrary to the conduction band Fermi sur- 
face, the valence band Fermi surface retains its spheri- 
cal topology for all t' /t ( see Fig. 1(b)). The experimen- 
tal Fermi surface of LaNiOa observed by Eguchi at al 
strongly resembles Fig.l(d).^^ 

The presence of large flat regions leads to nesting, and 
a tendency for CDW and/or SDW instabilities.^^ A sim- 
ple understanding of the effect of nesting is obtained from 
the Random Phase Approximation (RPA), in which the 
effect of interactions on the spin susceptibility is approx- 
imated by 
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where xq{^->^) is the non-interacting spin susceptibil- 
ity, and we took for simplicity a spin and momentum- 
independent interaction U . An instability is signalled 
by a divergence of x(0, k), which occurs on increasing U 
when the denominator in Eq. (2) vanishes. This occurs 
for the k which maximizes Xo(k) = xo(07k), which de- 
termines the wavevector of the spin ordering. In the case 
of perfect nesting, E^^ — £^q+k for every q on Fermi sur- 
face with the nesting vector k, and the non-interacting 
susceptibility is itself divergent at this nesting wavevec- 
tor, indicating an instability for arbitrarily small U . Al- 
though this is not true in general due to imperfect nest- 
ing, the flatness of the Fermi surface greatly strengths 
the tendency to instability. 

To check this directly, we calculate the zero frequency 
spin susceptibility, which in general in the Matsubara 
formulation is given by 



Xo(*w„,k) = (^kS'lk) 
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(3) 
(4) 



where the free electron Green's function is defined as 



Go{iuJn,q) = (4Cq) = {iuJn - Eq)- 



More details 

are explained in Appendix.A. Fig. 2 shows the calcu- 
lated zero frequency spin susceptibility Xo(k) as a func- 
tion of k for different ratios of t' /t. As expected, the 
spin susceptibility is sharply peaked at a particular cer- 
tain wave vector in the physical range of t' /t. Specif- 
ically, Xo(k) for t' /t = 0.15 shows the highest peak at 
k = Qn = 27r(|, |, |) which defines the nesting vector. 
Note that this is precisely the magnetic ordering wavevec- 
tor (in the cubic convention) observed in the insulating 
low temperature phase of the nickelates. Estimating the 
instabilty from Eq. (2), we obtain Uc ~ l/xo(Qn) ^ 2 
(see Fig. 2). 



III. HARTREE-FOCK THEORY 

A. Restricted Hartree-Fock Method 

Having established the nesting wavevector, we proceed 
to a (restricted) Hartree-Fock treatment of the ordering 
and MIT. We include interactions in the two-band model 
via an on-site Coulomb term U and Hund's coupling Jh^ 
defined from H = J^tb + ^int , 



JT.„t = UY,n^-JHj2s^ 



where rii = ^^ 
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As 
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discussed earlier, what is important here, because of 
the nesting physics, is the interaction between states 
near the Fermi surface. As such, the U and Jh terms 
may be thought of as simply a convenient parametriza- 
tion of the spin-independent and spin-dependent parts of 
these interactions, rather than literally in terms of atomic 
Coulomb and Hund's rule terms. 

To treat the problem in Hartree-Fock, we define a vari- 
ational wavefunction as the ground state of a fiducial 
mean-field Hamiltonian, which has the form of a non- 
interacting two-band hopping model plus linear "poten- 
tials" arising from coupling to SDW and CDW order 
parameters. Experimental results predominantly favor 
collinear magnetic ordering, of the form 



{Si) (X hi = iRe[^e*^--^^] 



(6) 



with complex variable ?/^ = \il)\e^^ . Fig. 3 shows different 
spin configurations which depend on the phase of 0. For 
instance, 6> = corresponds to "site-centered" spin or- 
dering in which the spin pattern is "up-zero-down-zero" 
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FIG. 3: Spin configurations depending on the phase of V^, 
^, along X axis, (a) shows "site-centered" spin ordering for 
^ = 0, (b) for intermediate — 7r/8 and (c) is "bond-centered" 
ordering for — n/A 



moving along a cubic axis, while = 7r/4 gives "bond- 
centered" ordering, and an "up- up- down- down" pattern. 
In the intermediate regime < 6 < 7r/4, the order is 
"off-center" as shown in Fig. 3(b). 

As already discussed above and in Ref.18, a CDW 
order parameter will be induced with Qcdw = 2Q^ = 
7r(l, 1, 1) as observed in experiment. This charge ordering 
is commonly known as "rock-salt" ordering and implies 
the electron density at site i is represented as 



(ni)ocpi = (-ir+J''+^'$, 



(7) 



where ^ is an Ising-type order parameter for the charge 
ordering. 

The full mean-field Hamiltonian from which the 
Hartree-Fock variational ground state is constructed then 
takes the form 
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(8) 
(9) 



The local exchange field h^ and the charge ordering pi 
couple to the spin operator S^ and the electron number 
operator n^ respectively. Note that we allow additional 
freedom in the variational state by letting the hopping 
parameters renormalize. That is 



Hth = Hti,[t^i,t' ^i% 



(10) 



The restricted Hartree Fock calculation proceeds by 
finding the ground state of H^ar'- 



Fvarl^o) = ^0|*0) 



(11) 



with the constraint of quarter-filling, i.e. one electron 
per site, ^- rii = N^ where N is the number of sites. 
The Hartree-Fock ground state |^o) is then a function 
of four dimensionless parameters: t' /t, \tlj\/t,^/t and 0. 
For each set of parameters, we calculate the variational 
energy 



Env = (*o|i^|*o), 



(12) 



which is then minimized over the dimensionless parame- 
ters, for fixed physical parameters t^t' ^U and Jh- 

To find l^o) in practice, we work in the reduce Bril- 
louin zone (BZ) determined by the four site magnetic 
unit ceh. We thereby end up with instead of two bands 
8 magnetic ones, constructed from the different pieces of 
the original BZ folded into the magnetic one. 



^nacxy^j — ^acxy^ ~r ^V^sdwy? 



(13) 



with n = 0, 1,2,3 (for four magnetic sublattices), where 
a is for two e^^ orbit als and a is for spin t^. In this basis, 

^tb = ^^^a5(k + nQsdw)4aa(k)Cn6c.(k) (14) 

k n 

The prime on k sum means the sum over the reduced 
BZ. In the same way, the density wave Hamiltonian in k 
space is represented as, 

i^dw = ^^^c|,+i^^(k)Cnaa(k)+h.C. 



k n 
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(15) 



We then find the single-particle eigenstates by diagonal- 
izing the 8x8 matrix of the variational Hamiltonian 
Eq. (9), and construct |^o) by filling the states up to 
the Fermi energy, determined by the requirement of 1/4 
filling. It is then straightforward to express £^hf in terms 
of the single-particle states and occupation numbers, and 
perform the minimization procedure (see Appendix. B for 
more details). 



B. Hartree-Fock Phase Diagram 

1. Two SDW states 

The resulting Hartree-Fock phase diagram for a typical 
situation with t' /t={).lb is shown in Fig. 4. We observe a 
metallic regime at small U and Jh , and two main ordered 
phases with stronger interactions. For large Jh^ site- 
centered SDW ordering with ^ = occurs, concurrent 
with strong charge order, generating an insulating state. 
This is natural because the large Jh favors pairing of 
electrons into spin S = 1 moments, requiring neighboring 
empty sites. More mathematically, such Hartree-Fock 
states minimize the Hund's term. For large [/, the bond- 
centered SDW with = 7r/4 occurs instead. This is 
again natural because the U term prefers uniform charge 
density, and with = 7r/4 in the cubic system (which we 
discuss here) no CDW order occurs. 

2. semi-metallic B-SDW 

Somewhat surprisingly, the bond-centered SDW state 
remains semi-metallic even at relatively large U within 
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FIG. 4: MIT bulk phase diagram for the model with cubic 
symmetry as a function of JJ jt and Jnjt^ with t' jt — 0.15. 
Here JJ is the on-site Coulomb interaction, Jh is the Hund's 
coupling and t is the nearest-neighbor hopping magnitude. 
The main phases which appear are a paramagnetic metal- 
lic state (M), metallic SDW (dark pink region close to M), 
insulating site-centered SDW (S-SDW), and semi-metallic 
bond-centered SDW (B-SDW). In between the S-SDW and 
B-SDW state, one observes an off-center SDW phase with 
< ^ < 7r/4 (pink region). The blue colored shapes show the 
points for which the optical conductivity is plotted in Sec. VI 



the Hartree-Fock approximation. Indeed, examination 
shows that the density of states is almost linearly van- 
ishing approaching the Fermi energy in this region, with 
a small non-zero value at E^^ which decreases with in- 
creasing U . This unusual behavior arises from the specific 
"up- up- down- down" magnetic ordering in this phase. To 
understand it, recall that the cubic lattice, viewed from 
the [111] direction, forms stacks of triangular lattice lay- 
ers. In the limit of strong bond-centered ordering, the 
spins on each triangular plane are fully polarized. More- 
over, electrons of one spin polarization are confined to a 
pair of parallel [111] planes, which together forms a hon- 
eycomb lattice when connected by the dominant nearest- 
neighbor hopping t. Thus in the limit of large U jt in the 
bond-centered SDW state, the appropriate tight-binding 
model is that of doubly degenerate e^ orbitals on honey- 
comb lattice. 



^0 = -E^#4aCi6 



(16) 



This model has four orbitals per unit cell due to the dou- 
bly degenerate e^ orbitals and the bipartite honeycomb 
lattice. Fig.5(a),(b),(d) and (e) shows the dispersion 
and the DOS of this tight-binding model for the cases 
t' jt = and 0.15. Without second- nearest-neighbor 
hopping t' jt = 0, the result contains two bands which 
are identical to those of the canonical nearest-neighbor 
tight-binding model for graphene, possessing two Dirac 
cones with linear dispersion at Fermi level. The similar- 
ity with graphene has led to the suggestion that such sys- 
tems might be used to engineer a topological insulator. ^^ 



With increasing t' /t^ the DOS saturates at a smah non- 
zero value approaching the Fermi level. This is because 
finite t' /t introduces both second- nearest-neighbor hop- 
ping and, more importantly, coupling between the honey- 
comb bilayers. The latter expands the Dirac points into 
small electron and hole pockets, in a similar manner as 
inter-layer coupling does in graphite. 



3. effects of orthorhombicity 

As discussed in Ref.18, the bond-centered ordering in 
the large U region is actually unstable to orthorhom- 
bicity (GdFeOa distortion) , which is present in all the 
nickelates save LaNiOa. This is expected on symmetry 
grounds to drive the SDW off-center. The off-centering in 
turn induces charge order. Thus at the symmetry level, 
when orthorhombicity is taken into account, the large U 
region is completely consistent with experiment. 

What of the metallicity in this region? In the 
graphene-like honeycomb bilayer, the Dirac point de- 
generacy is protected, as it is in graphene, by inversion 
symmetry. Inversion is indeed preserved by the bond- 
centered SDW in the ideal cubic system. It is, however, 
violated when both the SDW and orthorhombic distor- 
tion are present. Hence, we expect that orthorhombicity 
not only affects the centering of the SDW, it also tends to 
open a gap in the electronic density of states, converting 
the semi-metal to a true insulator. 

We now study this microscopically. A leading effect of 
the orthorhombic distortion is expected to be a crystal 
field splitting of the e^ orbitals at each Ni site. Therefore, 
we add the on-site orbital splitting term 



^ortho = ^ Di • cl^TabCib 



(17) 



Here we have suppressed the (diagonal) spin indices, and 
introduced Pauli matrices r in the orbital space. Us- 
ing the symmetries of the Pbnm space group of the or- 



thorhombic structure, we find (see Appendix. C) that the 
"orbital fields" D^ are all expressible in terms of a single 
vector D : 



D, 



((-i)^^+^^L)^, {-if^^y^Dy, D') 



(18) 



For simplicity, we consider this term in the effective hon- 
eycomb lattice model, Eq.l6, relevant for the large U 
case. Fig. 5 shows how the the DOS changes in the pres- 
ence of an orthorhombic distortion. A gap indeed opens 
for sufficiently large D, as plotted in Fig. 6. 



4. Limitations of the restricted HF theory 

Because we consider a restricted Hartree-Fock ansatz, 
some lower energy states that do not fit this ansatz 
may be missed in Fig. 4. For example, near the onset 
of SDW order, at relatively weak interactions, there is 
the possibility of an incominensurate SDW. This may 
be expected since the best nesting vector determined by 
the maximum of the susceptibility is not exactly at the 
commensurate value, but rather at Q ^ 0.47r(lll) (see 
Fig. 2). Generally, commensurate states are preferred 
at strong coupling, and if incommensurate phases ex- 
ist, they would be expected to change to the commen- 
surate ones with increasing interaction strength, via a 
commensurate-incommensurate transition. ^^ 

We have also neglected the possibility of spontaneous 
orbital ordering, which could occur in the cubic model at 
large U. Indeed, orbital degeneracy is crucial to the semi- 
metallicity found in the B-CDW phase, as we have seen 
above via the introduction of orthorhombicity. Sponta- 
neous orbital splittings (ordering) provide a mechanism 
for the cubic model to achieve a truly insulating state, 
which it must at sufficiently large U. However, we ar- 
gue that the absence of any observed orbital ordering or 
Jahn- Teller distortion is evidence that this physics is not 
relevant for the nickelates. 
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IV. STRONG COUPLING LIMIT 

The Hartree-Fock approach of the previous section is 
reasonable for weak to intermediate strength interac- 
tions, which we believe is most relevant for the more 
itinerant nickelates with R=Pr,Nd. For completeness, in 
this section we study the complementary limit of strong 
interactions, U/t^ Jn/t ^ 1. Here the two-band model 
is suspect, so the connection to experiment is less clear. 
However, we can at least qualitatively attempt to ad- 
dress the question of the interplay of charge and spin 
order in the strong coupling regime. Specifically, note 
that in the more insulating nickelates, with R=Eu,Ho,^^ 
charge ordering appears first upon lowering temperature 



from the paramagnetic metallic state, with magnetism 
occuring only at lower temperature. Thus it seems that 
in these materials there is a separation of scales, with 
the primary mechanism for the MIT being charge order- 
ing, and magnetism being secondary. In this section, we 
will see that this is indeed the case in one regime of the 
strong coupling limit of the two band model. The spe- 
cific parameters of this region do not, however, seem very 
physical, supporting the idea that in the more insulating 
nickelates a description beyond the two band model is 
needed. 

The strong coupling limit may be considered an ex- 
pansion in the hopping t^t' about the limit t = t' = 0. 
In the extreme limit, the behavior is determined entirely 
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FIG. 5: Dispersion (panels (a),(b) and (c)) and density of states (DOS, panels (d),(e),(f)) for the e^ tight-binding model on 
honeycomb lattice. In (a) and (d), t = 1 and t' = 0, we observe Dirac points with clear linear dispersion and corresponding 
linear DOS. The Dirac cone is stable to small t' — 0.15 as shown in (b) and (e). In (c) and (f), an orbital splitting induced by 
the orthorhombic distortion of the lattice is included, with t' = 0.15 and D = 1.5/a/3(1, 1,1). An induced gap is clearly seen. 
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FIG. 6: Plot of the single-particle gap A versus the orbital 
field |D|, for the honeycomb model with nearest-neighbor t — 
1 and t' — 0.15. Here we have arbitrarily taken the orbital 
field of the form D = (L>, D, D)/V3. 



In the remainder of this section, we will focus on this 
latter regime. Physically, we may consider the paired 
sites as bosons with spin S = 1. By introducing hop- 
ping perturbatively, we may introduce hopping and in- 
teractions between the bosons. In the perturbative treat- 
ment, we will, in addition to t/U^t/Jn <C 1, further as- 
sume tyt <C 1, which simplifies the algebra considerably. 
Below, we argue that the leading effects of hopping, at 
O(t^), induce charge ordering of the bosons, reducing the 
problem to an effective spin S* = 1 model. The spin de- 
generacy of the bosons is split only at the next non-trivial 
order, O(t^). This qualitatively agrees with the separa- 
tion of scales observed between charge and spin order in 
the nickelates. 



by the "atomic" Hamiltonian Hint in Eq. (5), which can 
be solved independently at each site, subject to the con- 
straint of proper total electron occupation (quarter fill- 
ing). There are two regimes, determined by the parame- 
ter a = U/Jh- For a > 1/4, the atomic ground state is 
one with one electron per site. In this regime every site 
is equivalent, and has four states available to it, due to 
the spin and orbital degeneracy. Further perturbation in 
t, f will therefore result in a spin-orbital Hamiltonian of 
the Kugel-Khomskii type. 

The other regime occurs when a < 1/4, and in this case 
the electrons prefer to segregate into two sets of sites with 
equal numbers in each: doubly occupied sites with total 
spin S* = 1, and empty sites. The ground state energy 
in this regime is Eq = —(1 — 2a)NJH^ where N is the 
number of sites. Here there are two sorts of degeneracies. 
First, foTt = t' = the location of the paired sites is un- 
determined, so there is a degeneracy of N\/[{N/ 2)1]'^ as- 
sociated with the different possible location of the pairs. 
In addition, for each of the paired sites, there are 3 spin 
states available. 



A. 0{t^ / Jh)'' charge ordering 

We first consider the effective Hamiltonian for the sys- 
tem at the leading non-vanishing order in perturbation 
theory, which is second order in hopping, for the case of 
a < 1/4. To formulate the perturbation theory , we treat 
the Hund's and Coulomb part as the unperturbed Hamil- 
tonian, Ho = i^int, and the hopping as the perturbation, 
1-Li = Hq. We denote the projection operator onto the 
ground state manifold of Ho at quarter filling by 7^. If 
1^) is an exact eigenfunction of the system with energy 
E^ then its projection into the ground state subspace, 
l^o) =7^1^) satisfies 
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where IZ = {Hq — E)~^ is the resolvent and Q = 1 — V 
. Eq. (19) is an implicit non-linear eigenvalue problem 
and we will only evaluate it perturbatively in Hi, then it 



becomes 
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« P7^i^(-l)"(7eQ7^in*o), (20) 



Thus, regardless of the specific spin states of the boson 
pairs, their nearest-neighbor interaction is always repul- 
sive for < a < |. The lower bound and hence charge 
order in the ground state follows. 



where to this order of accuracy, we can safely approxi- 
mate n^ {l-Lo-Eo)-^. 

The second order term in degenerate perturbation the- 
ory corresponds to n = 1 in Eq. (20), in which electrons 
make two consecutive virtual hopping transitions. The 
terms for three different types of hops can be combined 
(see Appendix. D for more details), up to an additive con- 
stant, into 
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Sj-l)]NiNj. (21) 



Eq. (21) gives the effective Hamiltonian at leading or- 
der for a < 1/4. To solve it, we note that Ni commutes 
with Hgff and is thus a good quantum number at ev- 
ery site. We then can easily see that the charge ordered 
states with A^^ = 0, 2 on the two rock salt fee sublattices 

0. 
Sj are 



of nil- 
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saturate a lower bound on the energy. 
This follows because, since the eigenvalues of 5^ 
bounded by -2, hence the effective boson-boson repulsion 
(the term in the square brackets in Eq. (21)) obeys 
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B. Magnetic interactions 



Notably, although the effective interaction V^^ ~ 
t^ / Jh determines the charge order in the ground state 
(and defines the energy scale separating it from uniform 
states), the spin degrees of freedom on the doubly oc- 
cupied sites remain undetermined at leading order. The 
spin physics is dictated by subdominant terms. Thus 
the appearance of charge order at a higher temperature 
than magnetism is a feature of this limit of the two-band 
model. Let us now consider the magnetic interactions in 
more detail. 

First we focus on the spin exchange between nearest- 
neighbor sites on the fee sublattice, i.e. second nearest- 
neighbor sites on the original cubic lattice. There are 
three lowest orders that we will consider: 0(t^/J^), 
0(tH'/Jl), 0{t'^/JH)- Although the effects of the t' 
hopping is a relatively small correction to the dominant 
t hopping, in the strong J limit, it is not negligible be- 
cause it can contribute at second and third order to the 
exchange between spins. Formally all these terms are 
on an equal footing if we take t' ~ t^ / Jh- We com- 
bine the contributions from different orders together (see 
Appendix. D for more details). The total spin exchange 
between nearest-neighbor sites on the fee sublattice is 
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Next we focus on the spin exchange between next 
nearest-neighbor sites. Calculation then shows 
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From the expression of Ji and J2, we obtain that if 
both t' and a are resonably small, there's ferromagnetic 
interactions between nearest-neighbors and antiferro- 
magnetic interactions between second nearest-neighbors 
on the fee lattice. Case (1) for 0{t^/Jfj) term is the 
dominant term for the ferromagnetic interaction. The 
negative sign for that case can be understood as arising 
due to the Hund's rule coupling on the intermediate k 
site, which prefers the two transferred virtual electrons 



to be in a triplet state. For the J2 exchange, however, 
because i^j^k are all along a single cubic axis, only one 
orbital can hop. For this reason, in the first hopping 
procedure (which is dominant) that contribute to J2 it 
is impossible to obtain a triplet intermediate state, since 
two electrons in a single orbital must form an antisym- 
metric singlet. This explains the antiferromagnetic sign 
of this exchange. 



Let us see what magnetic structure is expected from 
this exchange Hamiltonian. Since the fee lattice is a Bra- 
vais lattice, we can use the Luttinger-Tisza method to 
find the classical ground states. We simply Fourier trans- 
form the exchange couplings to obtain the energy of spiral 
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FIG. 7: Phase diagram of the classical ground state as a 
function of two dimensionless parameters a = U/Jh and 
/3 = tV^/t^ 



states with wavevector k. One finds 
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where x^ = cos/c^. Since this energy is quadratic in the 
x^, we can consider it as a quadratic form. The eigen- 
values of the form are 8(Ji + J2) and 8J2 — 4Ji (the 
latter is twofold degenerate). It is therefore positive def- 
inite if J2 > Max(— Ji, Ji/2). When this is satisfied, 
the minimum energy states are those with x^ = 0, i.e. 
k^ = ±7r/2. These are exactly the magnetic states ob- 
served experimentally. The phase diagram in Fig. 7 shows 
the classical magnetic ground state for different values of 
a and /3 = fJu/t^- We note that if t' = (/3 = 0) and 
Of = 0, the ground state appears to be ferromagnetic. 
When t'{j3) is included, the ferromagnetic Ji interaction 
is decreased, and the antiferromagnetic state will be sta- 
bilized. For the region a < 0.1, the magnetic ground 
state is antiferromagnetic when f ~ i^ / Jh- It is remark- 
able that one can obtain in this way the same magnet- 
ically ordered state as found from the itinerant nesting 
picture. 



C. Comparison with weak coupling limit 

According to the perturbation theory of the large 
Hund's coupling, charge order first appears at 0{t^ / Jh) 
and then magnetic ordering occurs due to perturbation at 
0{t^/jfj),0{t''^/JH) and 0{tH'/J%) . Since the magen- 
tic ordering arises from a temperature scale smaller than 
the charge ordering phase, this agrees with the exper- 
imentally observed intermediate charge ordering phase 
without magnetism. On the other hand, in the weak 
coupling limit, the charge ordering is always slaved to 
the primary magnetic ordering. ^^ 
V. CONFINEMENT EFFECTS IN THIN FILMS 

The success of the Hartree-Fock theory in reasonably 
predicting the charge and spin ordering in the more itin- 



erant nickelates undergoing a MIT suggests that the ap- 
proach may also be profitably applied to films. Recently, 
various growth issues have been overcome leading to epi- 
taxial films of good quality on several substrates with 
layer by layer control. One may expect that the MIT and 
related charge and spin ordering can be strongly modi- 
fied in thin films, due to both distortions (dependent on 
details of the substrate and growth conditions), effects of 
changes in chemistry at interfaces, and to quantum con- 
finement effects. Because of the difficulty of controlling 
the former two effects (which in any case are better stud- 
ied by first principles methods), we focus here entirely 
on the latter, and consider in this section the simplest 
possible model of a finite thickness film. That is, we sim- 
ply take the bulk tight-binding Hamiltonian and apply it 
to a finite thickness slab consisting of L unit cells in the 
confined direction, with effectively "vacuum" outside the 
slab, i.e. open boundary conditions. Given the impor- 
tance of Fermi surface shape in determining the nesting 
properties, we expect that quantum confinement alone 
can significantly modify the MIT properties and the or- 
dering in the insulating state. 



A. Single Layer, L — 1 



First of all, we consider the extreme case of a single 
Ni02 layer, following the methods used for the bulk. 
Here and throughout this section, we will neglect the 
symmetry-lowering effects that must be present in such 
a two dimensional structure, and in particular any tetrag- 
onal orbital splitting which is likely to be the dominant 
effect of this type. With this proviso, the Fermi surface 
and nesting properties are shown in Fig. 9(a). The two di- 
mensional Fermi surfaces show large flat regions similar 
to the bulk case. The zero frequency spin susceptibil- 
ity, Xo^(^)7 is shown in Fig. 9(e) (see Appendix. A). It 
is sharply peaked at Qg^w = 27r(|, |). Repeating the 
Hartree-Fock calculations for this case, using this SDW 
vector, we obtain the phase diagram in Fig. 8(a). The 
results are quite similar to the bulk case, except that the 
bond-centered SDW is insulating in this case, as the hon- 
eycomb lattice structure does not arise for a single square 
lattice layer. Somewhat surprisingly, the location of the 
MIT {U/t)c ~ 2 at J/f/t = remains almost unchanged 
from the bulk case. Naively, one would expect a decrease 
in {U/t)c in 2d, because the bandwidth is reduce by con- 
finement. We attribute the lack of such a decrease to 
decreased nesting in the two dimensional case, as can be 
seen by comparing Fig. 2 and Fig. 9: the susceptibility has 
a higher peak in bulk than in the single layer. 
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FIG. 8: Panel (a) shows the zero temperature phase diagram for a single layer L = 1 with the nesting vector Qsdw — 27r(|, |). 
As in the bulk case, the paramagnetic metallic phase (M) is stable for weak interactions. The dark pink region indicates a 
metallic SDW state, and the light pink region indicates an insulating off-center SDW. Panel (b) shows the phase diagram for 
three layers, L = 3, with the nesting vector Qfd^"^ = 27r(|, 0, 0). This nesting vector leads to a metallic B-SDW phase, which 
persists even for large U/t. The yellow region between the S-SDW and B-SDW phases is a metallic off-center SDW. 



B. Intermediate thickness films 

We now consider the intermediate cases with L > 2 
Ni02 layers along the z direction. In this case, the 
single-particle states can be taken as standing waves in 
the vertical (i) direction, with k^ = 7rl/{L + 1) where 
/ = 1,2,---,L. One obtains correspondingly 2L sub- 
bands (2 arising from the orbital degeneracy), each of 
which may have a Fermi surface. The calculated non- 
interacting Fermi surface and spin susceptibility for sev- 
eral values of L are shown in Fig. 9 (see Appendix. A for 
more details of the calculation of the spin susceptibility). 

From Fig. 9(g), we see that the peak of the susceptibil- 
ity varies considerably and in a non-monotonic fashion 
with L. While the case L = 2 (orange points in Fig. 9(f)) 
is quite similar to the result for the single layer, L = 3, 4 
are considerably distinct. For larger L, there is a slower 
variation of behavior, and by increasing the thickness 
to L = 30 (purple points in Fig. 9(f)), the bulk behav- 
ior (black line in Fig. 9(f)) is almost perfectly recovered. 
Thus we expect particularly distinct phase diagrams for 
the cases L = 3, 4, and focus on these below. 



The full Hartree-Fock phase diagram, assuming this 
wavevector, is shown in Fig. 8(b). Details of the calcu- 
lations for finite L, which are somewhat complicated by 
the many subbands, are given in App.B. Once again both 
site-centered and bond-centered SDW states appear, but 
the site-centered SDW occurs here only at very large val- 
ues of the Hund's coupling, Jn/t ^ 10, making it prob- 
ably entirely unphysical. Another distinction from the 
cases discussed previous is that the bond-centered SDW 
for L = 3 appears to be fully metallic. This is because the 
SDW with wavevector Q = 7r/2(100) describes stripes of 
electrons with all spins parallel in vertical stripes along 
the y direction. Thus the electrons are free to hop in 
this direction - actually they form "ladders" of two par- 
allel spin-aligned chains - and one has a sort of quasi- 
one-dimensional metallic state. Instabilities of the one- 
dimensional ladders would probably be expected beyond 
the Hartree-Fock approximation, and could lead to fur- 
ther charge/spin/orbital ordering and insulating behav- 
ior, but this is not within the scope of our study. 



2. L = 4 



L L = 3 

For L = 3, one observes comparable peaks in the 
susceptibility at two wavevectors: Q = 7r/2(100) and 
Q = 7r/2(110). The former is quite distinct from the 
ordering in the single layer and bulk cases. To decide 
amongst the two possibilities, we compared the varia- 
tional energy in the Hartree-Fock approximation for the 
two choices, and found that, over the full range of U and 
J, the total energy is lower for Q = 7r/2(100). Thus the 
model predicts quite distinct ordering in the trilayer case. 



One more noticeable feature in the spin susceptibil- 
ity plotted Fig. 9(f), is the large Q ~ peak for L = 4 
(see blue points). The Q = (uniform) susceptibility 
is simply proportional to the density of states, which is 
apparently enhanced for this film thickness. The origin 
of this enhancement is seen by inspecting separately the 
Fermi surfaces associated with individual sub-bands with 
discretized kz = 7rlz/{L + 1), shown in Fig. 10. One sees 
that the L = 4 case is unique in having three distinct 
Fermi surfaces (two hole and one electron) for the /^ = 1 
sub-band. Since the density of states is proportional to 
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FIG. 9: The zero frequency spin susceptibility, Xo(k), for finite thickness slabs of L layers, for the free electron tight-binding 
Hamiltonian with t' /t — 0.15. Plot (a) shows the purely two-dimensional single layer case. Here Xo(k) is sharply peaked at 
k = Qsdw — vr/2(ll). Plot (b) shows several cases with varying thickness with 2 < L < 30, compared with the bulk case 
L — oo. One sees that large L = 30 (purple dots) agrees well with the bulk susceptibility (black solid line). For smaller L, we 
see that the nesting properties change considerably. This is especially pronounced for L = 3 (green points), for which Xo(k) is 
sharply peaked at k = Q ?^ 7r/2(100), and for L = 4 (blue points), for which it is peaked at k ?^ 0. 
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FIG. 10: The Fermi surfaces for different subbands in the kx-ky plane with discretized kz = tiIz/{L + 1) for the cases L = 3,4 
and 5. The blue (red) lines correspond to conduction (valence) subbands. In the case L = 4, we see two valence subband Fermi 
surfaces at kz = 0.27r (Iz = 1), which is responsible for an enhancement of the DOS at Fermi energy. 



the Fermi surface area, this explains the observed en- 
hancement. Some understanding of this is obtained by 
inspecting the bulk Fermi surface, Fig. 1(d). It contains 
a large hole-like surface which has rather flat faces, par- 
allel to [001] planes. For the specific case L = 4 and 
L = 1, the discretized kz = 0.27r cuts across this rather 



flat when /^ = 1. As a result, the flat face, leading to the 
multiple two-dimensional sub-band Fermi surfaces. This 
enhanced density of states could potentially lead to fer- 
romagnetism for this case, but since ferromagnetism is 
notoriously over-estimated by the Hartree-Fock approxi- 
mation, we do not pursue this further here. 
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VI. DISCUSSION 

In the prior sections, we have studied a minimal two 
band model for the perovskite nickelates, with a focus on 



the MIT and the spin and charge ordering in the insulat- 
ing state. 
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A. Do Ave need the oxygen orbitals? 

In the minimal model used in this paper, we have elim- 
inated the oxygen orbitals to obtain an effective two or- 
bital Hubbard model. Several papers in the literature, 
however, claim that the oxygen states are crucial for the 
physics of the nickelates. Here we will discuss this issue, 
and argue that the importance of explicit inclusion of the 
oxygen states depends upon the questions being asked. 

In general, in the Fermi liquid paradigm, which applies 
to weakly to moderately correlated itinerant systems, the 
behavior of the electrons is dictated by the vicinity of 
the Fermi surface (s) only, and by the effective interac- 
tions amongst these states near the Fermi surface. The 
great insight of Landau in developing Fermi liquid the- 
ory was that the actual wavefunctions of these "quasi- 
particle" states are largely unimportant. Thus when 
it applies, any model which properly mimics the band 
dispersion near the Fermi surface (and its symmetry), 
and which captures sufficiently the interactions amongst 
the near- Fermi surface states, serves to correctly model 
the electronic behavior. It is well established now that 
LaNiOa, the metallic end- member of the RNiOa series, 
has a Fermi surface which is obtained from the inter- 
section of just two bands with the Fermi energy. These 
bands have Cg character, which can be mimicked by the 
minimal tight-binding model used in this paper. Pro- 
vided a band picture of the important electronic states 
near E^ is adequate, this basis is sufficient to describe the 
nickelates. The extent of the microscopic oxygen versus 
nickel character of the states is subsumed into the Bloch 
wavefunctions, which do not appear in the band Hamilto- 
nian^ and to a lesser extent in the effective interactions. 
We conclude that for low to intermediate energy prop- 
erties for which the two-band description is adequate, 
explicit treatment of the oxygen states is not important. 

However, one may ask questions - and conduct exper- 
iments - for which the oxygen states are obviously essen- 
tial. For instance, inelastic x-ray scattering can measure 
the relative fraction of Ni^+ and Ni^+ occupation of the 
Ni 3d states. Estimates for NdNiOs is that there is as 
much as 40% Ni^+. By neutrality, the Ni^+ can only 
arise through the presence of holes on the oxygen states. 
This implies the Bloch wavefunction associated with the 
the "oxygen bands" and "nickel bands" have in fact con- 
siderably mixed character. However, this does not affect 
the reliability of the two band model for the states near 
the Fermi energy. Indeed, the measurement of the Ni va- 
lence state is actually a measure of the occupied states, 
and hence is really related to the character of the filled 
valence band Bloch wavefunctions, not that of the near- 
Fermi surface states. Of course, by orthogonality, if the 
nominally oxygen states have mixed character, so too 
must the nickel states. 

Other high energy questions may be sensitive to the 
oxygen character. For instance, let us consider the prop- 
erties of an interface. In standard semiconductor sys- 
tems, an interface can be understood through band di- 



agrams, which include only the energies of the bands, 
and not their wavefunctions. Thus when this approach 
applies, the oxygen character is not important. In fact, 
band diagrams rely upon a semiclassical treatment which 
assumes that the electrostatic potential, carrier density, 
etc. vary slowly with respect to the lattice spacing. This 
in turn is correct in semiconductors due to their small 
effective mass and large dielectric constant. There is no 
need for this to apply to nickelate interfaces. 

In fact, it would be natural to expect a change in the 
oxygen character at an interface.^^ Consider an interface 
with a band insulator such as LaAlOs (LAO), in which 
there are no 3d orbitals near the Fermi energy. By neu- 
trality, in LAO the oxygen valence should be "exactly" 
(or at least much more so than in the nickelates) 0^~. 
This implies that the Ni 3d orbitals in the plane adjacent 
to the LAO are less able to hybridize with the intervening 
oxygens, since these states are "blocked". One can con- 
sider a simple model in which this physics is accounted 
for by ascribing an oxygen orbital energy e^/ for the inter- 
vening oxygens which is lower (so that here the electrons 
are more strongly bound to their oxygen) than the energy 
Cp for the same orbitals inside the nickelate, i.e. e^/ < Cp. 
The larger energy separation e^ — ep> > e^ — €p for the 
interfacial states implies reduces mixing of the nickel and 
oxygen states. Thus we expect that the Ni'^'^ character of 
the interfacial nickel ions should be reduced. As already 
remarked, this is a high energy property, related to the 
occupied states. However, the reduced mixing has im- 
plications at low energy as well. It implies reduced level 
repulsion between the 3d (specifically the d^^) and 2p 
states, so that the partially filled orbitals corresponding 
to the near Fermi energy states should be lowered relative 
to bulk nickelates near the interface. That is, the con- 
duction electrons feel an attraction to the d22 orbitals in 
the interfacial Ni02 plane. Note that, although oxygen 
physics induces corrections to its Hamiltonian parame- 
ters, the two-orbital model remains valid even for the 
interface. 

This physics may be relevant to recent experiments on 
LNO heterostructures. Several experiments have indi- 
cated the formation of an insulating state for very thin 
LNO films with only a few unit thickness. This appears 
at odds with the calculations in Sec. V, which find that 
the metal-insulator transition point is largely unchanged 
by confinement, even for very thin films. This model, 
however, neglects the induced orbital potential at the in- 
terface. One would expect this orbital potential to par- 
tially polarize the orbitals at the interface in favor of the 
dz2 states which conduct poorly in the xy plane. More- 
over, the shift of these orbitals renders inter-layer tunnel- 
ing non-resonant, which will further reduce the kinetic 
energy. Thus it is natural to expect the insulating state 
to be enhanced by this effect. In the future, we plan to 
investigate this in more detail by including the interfacial 
orbital attraction explicitly in the Hartree-Fock calcula- 
tion. 
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B. Strong versus intermediate correlation 



more insulating nickelates. 



In this paper, we have contrasted the limits of weak to 
intermediate correlation (and Hartree-Fock theory) and 
strong correlation (the perturbative approach in Sec. IV). 
It was argued that the strong coupling limit seems not 
very realistic. However, there are indications that some- 
thing beyond the weak coupling view is needed, at least 
for the more insulating nickelates, with R=Lu,Ho,Y. In 
these materials, the charge ordering and insulating tran- 
sition occurs above 500K but magnetism only sets in 
around lOOK. A factor of 5 or more discrepancy between 
these two scales is hard to reconcile with a weak-coupling 
picture. One type of strong-coupling picture is discussed 
by Anisimov et ai^^, in which the nickel charge state is 
regarded as Ni^+, which forms an 5* = 1 spin, while the 
mobile charge is actually in the form of holes on the O 
sites. The corresponding model would be a type of under- 
screened Kondo lattice. Charge ordering of the type seen 
in experiment is certainly possible, and would be viewed 
as the formation of collective Kondo singlets between 
two holes and a Ni^+ spin on half the lattice sites. ^^ To 
our knowledge, whether this actually occurs for a Kondo 
model of this type has not been established theoretically. 
This is an interesting problem for future study. A likely 
issue with such a Kondo description is that the band 
structure appears very different from the bands with Cg 
character predicted and observed in LaNiOa. Instead, 
the itinerant carriers must arise from oxygen bands, and 
it is not clear why this should in any way mimic the Cg 
structure. But perhaps the bands in LuNiOs etc. are rad- 
ically different from those in LaNiOs. If so, this should 
be testable experimentally. 

Some sort of intermediate coupling picture is also pos- 
sible. Indeed, even if the most insulating materials are 
at strong coupling, and, as we have suggested, PrNiOs 
and NdNiOa are better thought of in the SDW (weak to 
intermediate coupling) limit, then there are compounds 
in between. Here presumably a full description with all 
the orbital involved and charge fluctuations allowed in 
all orbitals is needed, and there is little simplicity to be 
found. Probably an approach which combines elements 
of ab initio theory and reasonable but ad-hoc treatment 
of interaction physics such as DMFT is the best in this 
regime. ^^ In this situation, it will unfortunately probably 
be difficult to identify any single mechanism for charge 
ordering. 

In our opinion, it is likely that one physical effect we 
have not so far discussed, the coupling to lattice phonons, 
is important. The Kondo singlet formation mentioned 
above would obviously benefit from a contraction of the 
neighboring oxygens around the Ni^+ spin in question. 
Indeed, it is this contraction which is actually observed 
experimentally, rather than any real electric charge den- 
sity. The same local phonon mode which would couple to 
the Kondo singlet would also favor charge ordering in the 
intermediate coupling view. It may be that this electron- 
phonon interaction gives a reasonable mechanism for the 



C. Experimental signatures 

It is desirable to understand how the different scenar- 
ios might be distinguished experimentally. We will focus 
here primarily on the expected consequences in the itiner- 
ant regime, as the primary focus of this work. However, 
we briefly discuss expectations for the strong coupling 
limits. In the strong coupling pictures, we would pre- 
sumably expect the insulating states to have a full gap 
to electron and hole quasiparticles. Moreover, local S = 1 
moments would be well-formed on half the Ni sites (form- 
ing an fee sublattice), prior to ordering into an antiferro- 
magnetic ground state. With these site-center local spins, 
it seems difficult to imagine an antiferromagnetic state 
with the symmetry of the bond-centered or off-center 
SDW, and we would expect a site-centered SDW (an- 
tiferromagnetic) order. This particular symmetry could 
be distinguished by a careful determination of local mo- 
ments at all the nickel sites from neutron or NMR//iSR 
measurements. 

Turning now to the itinerant regime, we consider the 
experimental consequences of the nesting scenario. First 
we discuss the thermal phase transition. In this limit, 
since the SDW drives the charge order, the two types 
of order should set in simultaneously at a single critical 
temperature. In Ref.18, it was shown that this transi- 
tion is theoretically expected to be first order for several 
reasons. These two observations are consistent with ex- 
periment. 

More detailed comparison can be made with electronic 
structure. We discuss in particular the implications of 
the nesting scenario for dc transport and optical mea- 
surements in the following. 



1. Transport anisotropy 

Transport is an important probe of the electronic struc- 
ture. In the nesting picture, the SDW order is di- 
rectly and strongly coupled to the quasiparticles, and 
hence should strongly influence the transport. The most 
qualitative feature of this coupling is that the SDW or- 
der imposes its lower lattice symmetry, and in particu- 
lar, spatial anisotropy, upon the quasiparticles. In con- 
trast, within the strong coupling view, the charge or- 
dering is dominant, and this charge ordering itself is 
not anisotropic (it doubles the unit cell but is compat- 
ible with cubic symmetry). We therefore expect that, 
when the nesting picture is valid, prominent transport 
anisotropy should be observed to set in for T < Tmit- 

We focus first on the bulk case, for which the SDW 
wavevector Qsdw = 27r(l/4, 1/4, 1/4) obviously breaks 
cubic symmetry. As discussed in Sec. IIIB2, the elec- 
tronic structure in the B-SDW phase is describable as 
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a set of weakly coupled honeycomb [111] bilayers, lead- 
ing (neglecting orthorhombicity) to a semi- metallic state. 
Hence we expect the B-SDW ordering to be accompanied 
by strong electrical anisotropy, with much larger conduc- 
tivity within the [111] plane than normal to it. 

We have calculated this conductivity at zero temper- 
ature using the Hartree-Fock quasiparticle Hamiltonian. 
From the Boltzmann equation within relaxation time ap- 
proximation, one has 



cr^^ = X^eVy 



8^^ 



-f'{en)]VnAk>nAk). (26) 



where r is a constant relaxation time, / is the Fermi dis- 
tribution /(e) = l/(e^(^-^) + 1) and Vn,^, = den{k)/dk^, 
where n is a band index. We have a total of 8 bands (2 
e^ orbitals x 4 magnetic sublattices = 8), and the band 
energies and velocities must be found numerically. Using 
k • p perturbation theory,^^ one has: 



^n,/2 



(k) = {i>„k\^^\^nk) 



(27) 



where H(/c) is the 8x8 matrix Bloch Hamiltonian. From 
the above formulae, we calculated the conductivity (j\\ 
parallel to the [111] axis and (j± normal to it. The ratio 
is plotted in Fig. 11 for the B-SDW state. As expected, a 
large anisotropy is observed once a significant magnetic 
order develops. 
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Electrical conductivity for L 
l,f = 0.15,6' = 7r/4 and $ 



3 and bulk with 
0. The isotropic 

Bulk 



conductivity in the bulk case (black line) , a ^ , decreases to 
zero with as a gap in the DOS develops with increasing SDW 
order. For the three layer case, L = 3, the conductivity shows 
a large anisotropy in the x — y plane once the SDW develops. 



effect due to SDW ordering should be observable below 

^MIT- 

For thin films, confinement effects may contribute to 
or modify the anisotropy. For instance, in the three layer 
case, we observed a change in the nesting wavevector to 
Ql=3 = 27r(l/2,0,0). In this state, the anisotropy axes 
imposed by the SDW are different. In particular, an "up- 
up-down-down" magnetic configuration along the x axis 
is stabilized, so that the spin polarized electrons are free 
to hop along y direction. Hence, in this case the low and 
high conductivity axes are the [10] and [01] axes, respec- 
tively. This is shown in Fig. 12, in which the magnitude of 
SDW, |?A|, is varied while fixing i = 1, ? = 0.15, = 7r/4 
and ^ = 0. Indeed, in this case the anistropic behavior 
is even more pronounced, for in the model the "hard" 
axis conductivity axx actually vanishes at T = in the 
limit of large SDW gap, while ayy saturates to a constant 
for arbitrarily large |?/^|, because the spin polarized elec- 
trons are free to hop along y direction. In this case, the 
formation of the SDW opens the Fermi surface. 



FIG. 11: Bulk conductivity anisotropy in the B-SDW state, as 
a function of the amplitude IV^I of the SDW order parameter. 



2. Optical conductivity 



Note that the same result would be expected to ob- 
tain for a thick film, where the behavior is predominantly 
bulk-like. In this case, the measureable quantity is the ef- 
fective two-dimensional conductivity tensor for the plane 
of the layer, which is usually an [001] plane. By symme- 
try, we expect the principle axes of the 2d conductivity 
to be the [11] and [1-1] directions, with different conduc- 
tivities along each in the SDW state. Note that in prac- 
tice this is complicated by the effects of orthorhombicity, 
which already should induce transport anisotropy even in 
the metallic state. However, we expect that this intrin- 
sic anisotropy is probably mild, and that a pronounced 



Optical conductivity is an other important probe of 
electronic structure. For LaNiOs, which is metallic at all 
temperature, experiment shows a reduced Drude peak 
compared to band theory, ^^ which may be considered 
as evidence of moderately strong correlation. However, 
apart from this quantitative renormalization of the low 
energy Drude part, the theoretical optical conductivity 
obtained from the simple two eg band model reproduces 
experiment fairly well up to uj ^ 2eV }^ Applying the 
same analysis to the magnetically ordered phases in our 
bulk phase diagram. Fig. 4, we obtained strikingly differ- 
ent results as a consequence of SDW formation. 



The calculations are made using standard linear re- 
sponse theory within the Hartree-Fock variational Hamil- 
tonian. From the Kubo formula, the real part of optical 
conductivity (Ja(3{^^'k) is related to the imaginary part 
of current-current correlation Ua/si^^'k):^^ 



^al3 



(Pk 

(2^ 



noe 



n,^(i],k) + ^^fo/3 (28) 



m 
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with wave vector k, frequency Vt^ average density no and 
electron mass m. The current-current correlation func- 
tion with imaginary frequency iVti is defined as 



Iic.^{inu k) = ;^ ^ jf{\^)jf{\^) ^ Y. GadiiuJn + i^i. k)Gcb{iuJn. k). 



vol 



abed 



At zero temperature, this can be calculated from the spectral representation (see Appendix. E), 

Im[n,^(l],k)] = ^C*^m'C'^m/^^m(c^)A^'(c^ + ^)(nF(c^)-n^(c^ + l])), 



(29) 



(30) 



r 



where A^(cc;) = j/Kuj — E^n + M^)^ + 7^], with 7 a 
small scattering rate (imaginary part of the first order 
self-energy correction Im[I](cJn)] = —i^sgn{ujn)) added 
by hand, and 0J^(k) is the a component of TTith eigen- 
state. riFiuj) = l/(e^^ + 1) is Fermi distribution. 

Fig. 13 shows the optical conductivity calculated in this 
way for each of the different phases (taken at the spots 
marked by symbols in the phase diagram in Fig. 4). The 
above-mentioned comparison of theory and experiment 
for the paramagnetic metallic state is shown in panel 
(a), taken from Ref.ll. The development of SDW or- 
der strongly suppresses the Drude peak, as expected, 
which can already be seen in the metallic SDW state 
when the density of states at the Fermi energy is still 
non-zero (but small). Fig. 13 (b). Interestingly, a small 
peak appears instead at u/t ~ 0.3. This peak arises 
from a transfer of spectral weight from low frequency to 
above the SDW gap. Panel (c) shows a{uj) for the B- 



SDW state, which has a semi-metallic band structure. 
One observes a linear increase of Re[(j(co')] for small fre- 
quency cj, which is similar to the behavior expected from 
the Dirac points in graphene, and indeed arises from the 
honeycomb [111] bilayer structure of the spin-polarized 
regions, as discussed in Sec. IIIB2. Here we have plot- 
ted the powder- average conductivity, since the full ten- 
sor is anisotropic as discussed above. This calculations 
has neglected orthorhombicity, which would introduce a 
gap at low energy and thereby interupt at least part of 
the linear region. However, a linear increase of (j{uj) at 
low frequency was indeed seed in bulk experiments on 
NdNiOs below the transition temperature.^^ Finally in 
Fig. 13(d), we plot the optical conductivity for for large 
Hund's coupling J, in the S-SDW where strong charge 
order is present. A large gap opens in the spectrum, 
resulting in zero Re[cr(co')] up to uj/t ^ 1. 



r 



D. Summary 



We have presented a theoretical analysis of the metal- 
insulator transition in the nickelates from a minimal two- 
band model and Hartree-Fock theory, which we argued is 
appropriate for the itinerant limit of weak to intermedi- 
ate correlation. This picture of the metal-insulator tran- 
sition can be tested in various ways, as suggested above, 
and appears to us to be the most consistent one for the 
materials NdNiOs and PrNiOs, located close to the zero 
temperature MIT phase boundary. For the more insu- 
lating nickelates, a different type of theory is required, 
involving stronger correlation and possibly an important 
role for electron-lattice coupling. Both further theoretical 



work in clarifying the mechanism for the MIT transition 
in those materials, and experimental work which can test 
the itinerant picture (such as measurement of transport 
anisotropy), would be very desireable. Finally, we have 
shown that quantum confinement alone cannot explain a 
Mott insulating in ultrathin LaNiOs films, and suggested 
a physical mechanism by which the observed insulating 
state might obtain. It will be interesting to pursue this 
question further in the future. 
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FIG. 13: The real part of the optical conductivity, a{(jo), for each phase (colored star) in the bulk phase diagram (see Fig. 4). 
In (a), the paramagnetic metallic phase shows a large Drude peak and a small hump (see inset plot). The hump is related to a 
region of large DOS for interband transitions (see Ref.ll). In (b), the metallic SDW phase has a reduced but non-zero Drude 
peak and a small second peak due to the SDW gap. Plot (c) shows the case of the semi-metallic B-SDW phase, for which a 
linear increase of cr(uj) for small frequency u is found, related to the linear dispersion near the Fermi level. It also shows strong 
anisotropy between the conductivity crfm] (along the [111] direction) (blue line) and crfm]^ (perpendicular to [111]) (green 
line). In plot (d), a large gap is visible in the S-SDW phase. 
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Appendix A: Dynamical Spin Susceptibility for Free 
Electrons xo{^,k) 

In this section, we derive the dynamical spin suscep- 
tibility for free electrons both for bulk and finite layers. 
In general, dynamical spin susceptibility for Matsubara 
frequency iQn^ and wave vector k can be represented as 
following 
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Xo{i^n,k) 






N 



;^(X1 XI 4CraCr/3C;,„,Cr'/3'Cr;^CT;,^,e 



r"! A" r.r'' . 



^ .rr^ . ^*k(r-r') 



rr' q;/3,q;'/3'' 



{qi} rr' olol' 



N 



{qi} aa' 



■^(11 I]4i"Ck+qxaC^3a'C_k+q3a'(-l)"+"') 



4Ar 



qiqs aa' 



/da 1 
-^- ^ Tr[Go(iw„, g)Go(i(w„ + f^„), g + fc)] 



2 7 (27r)3;9 



(Al) 
(A2) 

(A3) 

(A4) 

(A5) 

(A6) 

(A7) 



First of all, k(qi) is four dimensional vector which in- 
cludes both Matsubara frequency iQni'f^^n) and the wave 
vector k{qi) in three spatial dimension. In the same 
way, r includes both imaginary time r and spatial di- 
rection r. Qf,/3,a' and f3^ are for spin t,| and a repre- 
sents Pauli matrix and Sr = Xlra/3 ^ra~i^^r/3 (with ig- 



noring orbital indices for simplicity). Fourier transform 

Sk = ^ Er Sre*^"" and cj^ = ^c^c^e^'^'^, free electron 

Green's function Go{iuJn^q) = (c^Cq) = {ioUn — Eq)~^ . 
From the last equation of Eq. (A7), we sum all the Mat- 
subara frequenies using the following trick 



(5 ^-^ iuJn — X i{uJn + ^n) " x' iftn ^ X — X 



-{nrix) -nrix)) 



(A8) 



where Fermion distribution is defined as npix) = 
l/(e^^ + 1). For simphcity, we represent doubly degener- 
ate eg orbitals tight-binding model using Pauli matrices 



(T^ Htb{k) = eo{k)l + e{k) • a. Then finally analytic con- 
tinuation leads 



Xo 



^ ' ^ 2 J (2^)3 Z^ n^xa{q)-x,{k^q) 



a,be± 



{np{xa{q)) - np{x,{k + q)))^{l + ab^^^^j^±^^) (A9) 



where x±{k) = eo{k) ± |e(/c)| — ja. Here, N± is the number of sites on its perpendicular x — y 

For finite layers (along z direction), it has discretized plane. The spin susceptibility Eq. (Al) is represented 
kz = ^lz/{Nz + 1) where U ^ {1, 2, 3 • • • N^} and N^ is 
the number of layers. 



N^N, 



Y^ciJ'^-^-sink. 



(AlO) 



18 



9^ {_^\oi-\-<^' , II II, 



(^^^-^) = i^(EEE 



-^(EEE 






smqlvz sin(/2?'« sin(/|r^ sinq|r^e '^""^ ''»-') 
(-1)"+"' t t 



(All) 



^iV,( ^ 5(agf + bq^ + fc,))iV,( ^ 5(agf + bql - k,))) (A12) 
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o6e± 



a6e± 



A/-. \ 2^ Z^ 2^ 2^ '^9iL9ja'^fe-L+9i^92"'^g^g|a''^-'=-L+93'94"' 



N^N, 



1^13 ilt} ""' 



( J2 H^lt + bq2 + kM J2 ^('^ll + bql - K))) 



abe± 



abe± 






{S(ql + ql + fc.) + 5{ql + 9| - k,) - S{ql - g| + fc.) - 5(9? - g| - fc,))^ 



(A13) 



(A14) 



From Eq. (All) to A13, we abbreviate Matsubara fre- 
quency indices for simple representation. 6 functions in 

Eq. (A14) can be rewritten S{qi + gf + ^z) = ^{h + 
h + lz){mod{2{N, + 1))) where q^ = 7rk/{2N, + 1) and 



Appendix B: Detailed Hartree-Fock Calculation 

Our variational Hamiltonian Eq. (14) can be diagonal- 
ized by writing 

C„aa{k) = J2^n:ik)CAaik), (Bl) 



with A = 1 . . .8 indexes the eigenstates for each k. With 
appropriate choice of ^ the diagonalized Hamiltonian be- 



^var = ^^e^(fc)cL(fc)c^Jfc). (B2) 



k Aa. 



Here we have used that e^(/c) are independent of a. 
This can be seen since the transformation Cnaa{k) -^ 
{ — l)^Cnaa{k) maps a -^ —a. Hence 



<i>^-{k) = i-ir^^:ik), 



(B3) 



and the energies are independent of a. Now we take the 
expectation values of each term. Ground state |^o) is 
nothing but occupying all the quasiparticles states be- 
low the Fermi energy. First, we consider the expectation 
value of J^tb • 



(^tb) = E E E E ^-^"(^ + "^) i^l^naik))* cl^^nk)np{eA{k)). 



(B4) 



k n ab A 



19 

where ni?(e) is the Fermi function. Next, consider the expectation value of on-site Coulomb interaction Hu = U J2i ^f • 

i i 

= N Yl 4a(^l)Caa(^2)4^(A:3)c5^(A:4)4i+/c3,fe2+/c4 

I ' 

= TV ^ ^ 4iaa(^l)Cn2aa(^2)c|,35^(A:3)c^^5^(A:4)4i+fc3,fe2+fc4^ni+n3,n2+n4(mod 4) 
{ki}{ni} 

= ^EE E (Ca(^l))*<^afe)(C(M)*C(M 
{ki}{ni}ABCD 

><cL(^l)cBa(^2)cJ.^(^3)c^;3(^4)4i+/c3,/c2+/c4^ni+n3,n2+n4(mod4). (B5) 



Now we take the expectation value. There are both 
Hartree and Fock terms. 



J 



E-0 = ^EEEEE[«a(fci))*ea(fci)(CM*C1(M 

\ i I /ci,/e3 {ni} A,C ah a(3 



Using Eq. (B3), Eq. (B6) can be simplified after summing 
the spin indices a, /3 



E^n = ^ E EEE[2(i + (-ir^"^)«a(fci))*<+(fci)(C6(^3))*c.(fc3) 

\ i / ^1^3 {rii} AC ab 

-2 «+ (fci))* <^^+(fc3) (C3t(fc3))* <t(^i)]"F(eA(fci))nF(ec(fc3))<5„,+„3,„.+„.(mod 4) (B7) 



In the same way, the expectation value of Hund's cou- 
pling is represented by 
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ESf) = ^EE E [(<a(A^l))V«/,0^f„(fc2)]-[(Cl-'(fe))*'^"'/3'<'l'(M]nF(eA(fcl)W(ec(fc3)) 
\ i I ki {ni}ABCD 



4^- 



kik3 {m} AC a(5a'(5' 
nF(eA(/^l))nF(ec(^3))^ni+n3,n2+n4(mod 4) 



4^- 

kiks {m} AC ab 
nF(eyl(fci))nF(ec(fc3))<5ni+n3,n2+n4(mod4) 

I 



(B8) 



Finite layers N^ along i direction lead the discretized 
k^ = itIz/{L + 1) where / = 1, 2, 3 • • • A/"^ 



/v EE^ 

k\kz n 



,{k±kz)e 



lk±ri_i 



smkzVz (B9) 



Now, the interaction term Eq. (B6) has modified func- 
tion of 5fc^+/e3,fc2+fc4 ^or both Eq. (B6) and Eq (B8) 
need l/4(^;,i^+;.^^;.^+;.^(4 + 2^2(fci+fei)("^od(27r)) + 
2(52(/ci-/c3)(mod(27r)) The last two S functions corre- 
spond to ^(^i+^3),L+i and ^^^,^3- 



Appendix C: Orthorhombic GdFeOs distortion 

In this section, we discuss the point symmetries of 
the orthorhombic lattice (GdFeOs type perovskite) and 
study how this symmetry operators constraint on-site 
splitting vectors D^ defined in Eq. (17). We first the de- 
fine four basis sites of the orthorhombic lattice (in cubic 
coordinates): 



ri 
r4 



(0,0,0), 
(1,0,0), 
(0,0,1), 
(1,0,1). 



(CI) 
(C2) 
(C3) 
(C4) 



The orthorhombic space group has three point group op- 
erations (in cubic coordinates) : 



Pi 

P2 

Ps 



{x,y,z) — > {l-y,-x,-z), (C5) 

{x,y,z) — > (-1 -X, 1 -?/, 1 + z), (C6) 
{x,y,z) — > {-l-x,l-y,-z). (C7) 



One finds that Pi interchanges sites ri o r2, and rs o 
r4, while P2 interchanges sites ri ^ rs, and r2 o r^. 
The inversion Ps leaves the basis unpermuted. Taking 



the usual cubic basis of dx2_y2 and ^^2 orbitals, one then 
readily finds the transformations of creation/annihilation 
operators: 



Cl 


-^ 


-r"C2 


C2 


-^ 


-r^ci 


cs 


-^ 


-r^C4 


C4 


-^ 


-r^cs 



Pi'^i. Lz: P2'^<:' ^. 7 ^3:x, 



ci -^ Cs 

C2 -^ C4 

Cs -^ Ci 

C4 -^ C2 

(C8) 
where the last equation indicates that Ps acts as the iden- 
tity in both the orbital and sublattice space. From this 
we see that Ps places no constraints whatsoever on the 
orbital fields. Invariance under the first and second trans- 
formations then allows all four orbital fields to be deter- 
mined from one. One finds: 



r>2 = {-D\-D\D^), 
Di = {-D'',-Dy,D^). 



(C9) 
(CIO) 
(Cll) 
(C12) 



Thus there are two and not four different orbital fields 
appearing. Taking into account the coordinates of these 
basis sites, we can finally write a simple form which is 
basis independent: 



D, 



{^{-ifi+yiD'', {-if'^y'Dy, D^ 



(C13) 



Appendix D: Degenerate perturbation theory 
calculation in the strong coupling limit 

1. 0{t^ / Jh)'* charge ordering 

There are three possible types of hops at second order: 

1. An electron hops from a double occupied site to an 
empty site, and then back. This lowers the energy 
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when occupied sites are adjacent to empty sites, 
and so results in an effective repulsion between bo- 
son pairs. 

2. Both electrons from a doubly occupied site hop 
onto the same, previously empty, site. This results 
in an effective hopping of the bosons. 

3. In the case where neighboring sites are occupied 
with bosons, there can be exchange if the spins of 
both bosons are not parallel. 

The terms in the effective Hamiltonian corresponding to 
the above three procedures can be written as 



0/(1) 

^eff 

(2) 
eff 

(3) 
eff 



^bj 






n 



ij iaa jba 
lb J 



,t 



(Dl) 
-2'Pn!cLc^kaT^Qt'^ic],pC,^p'P. (D3) 



rttcir^^^nQt':,fci,c,,,v 



^ij ^iaa^jba 



ji jba iaa ' 



ij^ic(3^jd(3' 



All three terms include implied sums over nearest- 
neighbor sites i and j. Here we have neglected 0{t''^) 
contributions which are parametrically small in the limit 
considered. The factor 2 in Eq.(23) arises from the fact 
that electrons can hop first from site i to j or vice- versa. 
Using the exact form of the hopping matrix in Eq. (1), 
one finds that the second effective Hamiltonian vanishes, 
^eff = 0. This can also be understood from simple or- 
bital considerations: only one of the two orbitals over- 
laps along any of the principle directions. Since both 
electrons must be transferred for the pair to transfer, the 
boson hopping vanishes. Due to the absence of the pair 
hopping, the effective Hamiltonian commutes with Ni. 



2. Magnetic interactions 

Consider a plaquette on the original cubic lattice, we 
name the occupied sites i and j and empty sites k and j, 
such that i and j are next nearest-neighbor on the square 
plaquette. We calculate terms for nearest neighbor spin 
exchange at different orders one by one as follows. 

1. o^tVJ^H)-- 

(a) One electron from each of sites i and j hops 
to site /c, and then the two electrons at site 
k return to i and j. There are four distinct 
time orders in which this process can occur 
and they contribute equally. The same proce- 
dure can also happen to sites i, j and I. This 
give the coefficient 8 in front the Hamiltonian 
below 



0/(4) 
^eff 



^^'^jk^jc-f^kf-f^Q^ik^iaX^keX 



nattc 



dc J\ 



kj kda jca 



C,. 



yf'^^ti-CuhoC-oJ . 



ki'-kb(3'-ia(3' 



(D4) 



(b) One electron at site i hops to site k and then to 
j , it forms a singlet state with another electron 



at site j, and then one of two electrons forming 
a singlet hops back to site k and then to i (For 
brevity, we will not write down the effective 
Hamiltonian of the other hopping procedure 
from now on) 

(c) One electron from site i{j) hops to site k{l), 
so now four corners of the plaquette are all oc- 
cupied with single electrons, then the electron 
at site k{l) hops back to j{i) 

(d) One electron at site i hops to site k and then to 
jf , it forms a singlet state with another electron 
at site j, and then one of two electrons forming 
a singlet hops to site / and then to i. 

Combining the four terms, we have 
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(D5) 



(a) One electron at site i hops to j via next 
nearest-neighboring hopping, it forms a sin- 
glet state with another electron at site j, then 
one of two electrons forming a singlet hops 
back to site k and then to i 

(b) One electron at site i hops to /c, another elec- 
tron at site j hops to site i via next nearest- 
neighbor hopping, then the electron at site k 
hops to j. 

Together these two terms give 



Jl 



(2) 



tH' 



J'h 



l-4a 
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5 + 4q; 1 - 4a 



(D6) 



3. Oit'yjH): 

One electron at site i hops to an fee nearest- 
neighbor j, forming a singlet state with another 
electron at site j, then one of the two electrons 
forming the singlet hops back to site i. We obtain 



Jl 



(3) 



f- 



JH5 + 4a 



(D7) 



The spin exchange coupling between nearest neighbor 
is then Ji = jf ^ + jf ^ + jf \ 

For second nearest-neighbor spin exchange, consider 
three sites i,k and j along the same cubic axis, where 
i{j) and k are nearest neighbors of the original cubic lat- 
tice. Sites i and j then correspond to the second nearest- 
neighbor sites. Then there are two possible ways of the 
hopping procedure, which is of identical hopping order 
to the first two cases of 0{t^/Jfj) terms. 
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Appendix E: Optical conductivity 

The current-current correlation function with imagi- 
nary time is defined as 



n„^(r,k) = -^{Tr-r'ji{r,k)j^{T',k)) = A ^ jab(k)j-crf(k)G„,(r' - r,k)Ge6(r - r',k) 



vol 



abed 



where j^^ = dH^^ /dka, Gad{^' — ^^^ is a retarded vector k and a prefactor 2 for spin sums. Fourier trans- 
Green's function with imaginary time r' — r and wave form with Matsubara frequency iujn leads 



n„;3(if},,k) = -^^j2''(k)i^''(k)^^Gad(iw„+iO;,k)Ge6(iw„,k) 



(El) 



abed 



Green's function Gadi'i^^n-,^)^ 
G„d(iw„,k) = J dTe{T)Gad{T,k)e'^"^ 



= E 

mm' 

= E 



-/3e„ 



-/3e, 



lUJ, 



-{m'\cd\m){m\cl\m') 



n \^ ^rn i^ 



6-(k)0i(k) 



The last term in Eq. (E2) is for zero temperature with 
the imaginary part of the first order self-energy correc- 
tion Im[i;(a;n)] = —i'ysgn{cjOn)^ the a component of m 
eigenstate 0^(k). By substituting Eq (E2) to Eq. (El) 
and using analytic continuation iQi -^ Q -\- irj, the imag- 
inary part of the current-current correlation function is 
represented as 



iujn + i'jsgn{ujn) - {Em - IJ^N) 



(E2) 



lm[Ua^{n,k)] 



E^-^-'C'^l 



duj 



Am{c^)Am'{^ + ^){nF{c^) - UpiuJ + ^)) 



(E3) 



where Am(oj) = ^/[{oo - Em + /nN)'^ + 7^] and Fermi 
distribution nF{uj) = l/(e^^ + 1). 
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